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Abstract 

We analytically work out the orbital effects caused by a Rindler- 
type extra-acceleration A^in which naturally arises in some recent mod- 
els of modified gravity at large distances. In particular, we focus on the 
perturbations induced by it on the two-body range p and range-rate 
p which are commonly used in satellite and planetary investigations 
as primary observable quantities. The constraints obtained for Aftin 
by comparing our calculations with the currently available range and 
range-rate residuals for some of the major bodies of the solar system, 
obtained without explicitly modeling ^Rin, are 1 — 2 x 10"^'^ m s~^ 
(Mercury and Venus), 1 x 10^-^'' m s~^ (Saturn), 1 x 10^^'"' m s~^ 
(Mars), while for a terrestrial Rindler acceleration we have an upper 
bound of 5 X 10~^^ m s~^ (Moon). The constraints inferred from the 
planets' range and range-rate residuals are confirmed also by the latest 
empirical determinations of the corrections Aw to the usual Newto- 
nian/Einsteinian secular precessions of the planetary longitudes of per- 
ihelia zu: moreover, the Earth yields < 7 x 10^^® m s^^. Another 
approach which could be followed consists of taking into account Amn 
in re-processing all the available data sets with accordingly modified 
dynamical models, and estimating a dedicated solve-for parameter ex- 
plicitly accounting for it. Anyway, such a method is time-consuming. 
A preliminary analysis likely performed in such a way by a different 
author yields ^ < 8 x 10^^"* m s^^ at Mars' distance and A < Ix 10^^** 
m s~^ at Saturn's distance. The method adopted here can be easily 
and straightforwardly extended to other long-range modified models 
of gravity as well. 

Keywords: Experimental tests of gravitational theories, Modified theories of 
gravity, Celestial mechanics, Ephemerides, almanacs, and calendars. Orbit 
determination and improvement 
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PACS numbers: 04.80.Cc, 04.50.Kd, 95.10.Ce, 95.10.Km, 95.10.Eg 



1 Introduction 

Recently, Grumiller [T] has constructed an effective model for gravity of a 
central object of mass M at large scales. Starting from the most general 
space-time metric endowed with spherical symmetry in four dimensions [2] 

{dsf = gAB{x'')dx'^dx^ + ^^{x'') [((i^)^ + sin^ 9 {d^f] ,A,B = 0,1, (1) 

where is a 2-dimensional dilaton field, it is possible to reduce the 

4-dimensional Einstein-Hilbert action to a 2-dimensional dilator^ one [3]. 
By writing down the most general 2-dimensional action for qab and ^ 
consistent with spherical symmetry and with additional assumptions like 
power-counting renormalizability, analyticity, etc., and by considering the 
long-distance limit, an additional term with respect to the usual ones of the 
General Theory of Relativity (GTR) appears. It is [l] 

gABix^dx^dx'' = K\dxy - ^ = r, (2) 

with 

K^ = l-'^-Ar^ + 2Ar. (3) 
r 

While M and A are parameters retaining their usual GTR meanings. A, ac- 
cording to its geometrical interpretation [3] , yields a radially directed Rindler 
extra-acceleration Aj^^ = Ac? (RIN in the following) [2j. 

Since in Ref. [1] it has been envisaged the possibility that RIN may 
explain the anomalous behavior exhibited by both the Pioneer spacecrafts 
after they passed the threshold of about 20 astronomical unit (au) [U [5], 
traditionally expressed in terms of an unexplained radial acceleration 

Apio = (8.74 ± 1.33) X 10"^° m s'^, (4) 

it is worthwhile to study in more details the phenomenology of RIN in 
the solar system. Incidentally, let us recall that recent works [61 [7] point 
towards an explanation of the Pioneer anomaly in terms of mundane, non- 
gravitational effects peculiar to the probes. On the other hand, if the Pioneer 



^Dilaton gravity models in two dimensions have recently attracted much attention 
because of their implications for a list of well-known problems in quantum gravity [3] . 
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anomaly was really a genuine dynamical effect of gravitational origin, it 
should also affect the orbital motions of the solar systems bodies moving in 
the space regions in which it manifested itself in its presently known form. 
For such an issue see, e.g., Ref. [8] and references therein. 

In Section [2] we analytically work out some effects caused by RIN on test 
particles' orbital motion around a central body. In particular, in Section 
12.11 we deal with the two-body range, while in Section 12.21 we treat its time 
derivative, i.e. the two-body range-rate. In Section [3] we put constraints on 
RIN by comparing our results with the most recently produced range and 
range-rate residuals for several solar system bodies. Section|3]is summarizes 
our findings. 

Finally, let us note that the strategy presented here can, in principle, 
be straightforwardly extended to other exotic effects predicted by different 
long-range modified models of gravity. 

2 Analytical calculation 

In principle, simple back-to-the-envelope computations may be performed 
by noting that, given a constant and uniform extra- acceleration ^Rin acting 
on a test particle in orbital motion around a central body, its position and 
velocity shifts caused by ^Rin are roughly given by the product of vlRin by 
the second and first powers, respectively, of a characteristic time T of the 
system considered which, in the present case, is the particle's orbital period 
Pb- Anyway, such a naive approach is not able to tell us if such a kind 
of perturbation does actually affect the orbital motions with non- vanishing, 
long-term effects. Moreover, also by a-priori supposing that it is just the 
case, the correct order of magnitude of them may not be correctly inferred 
because of the true details of the orbit like, e.g., its eccentricity e. Indeed, 
since it is usually very small for typical solar system bodies, its presence or 
its absence in the expression for a certain orbital effect may substantially 
alter its size. Thus, it is mandatory to explicitly work out in full details the 
perturbations induced by A-^^^ on some features of test particles' orbits: we 
will choose th^ two-body range p and range-rate p because they are direct, 
unambiguous and accurate observables very common in Earth-satellite and 
Sun-planets investigations. In Section [2. II the perturbation Ap is computed, 
while Section 12.21 is devoted to Ap. 

^We refer to two test particles A and B orbiting the same central body and reciprocally 
interconnected by means of some artificial, man-made electromagnetic links. 
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2.1 The two-body range 

The orbit of a test particle around a central mass in presence of a dynamical 
perturbation A of the usual two-body, pointlike Newtonian monopole can 
be written as 

fp = fu + Af. (5) 

The subscript!! P denotes the perturbed orbit, while U labels the unper- 
turbed Keplerian ellipse. To avoid possible misunderstandings, we remark 
that what will, actually, be compared to the data in Section [3] is only ob- 
tained from the perturbation Af. As a consequence, standard perturbative 
approaches like those relying upon the Gauss or the Lagrange equations for 
the variation of the Keplerian orbital elements [13] can be adopted to work 
out the orbital effects of A. They make use of the Newtonian trajectory 
as unperturbed, reference path. In principle, it is possible to assume as 
reference orbit a fully post-Newtonian one [9l|T0], and work out the effects 
of a given small extra-acceleration A with respect to it according to the 
perturbative scheme set up by the authors of Refs. [9l [TO]. It is a general 
relativistic generalization of the standard perturbative approach based on 
the planetary Lagrange equations [13]. In the present case, given the ex- 
tremely tight bounds on ^Rin which will be inferred in Section O we may, in 
principle, safely apply such relativistic perturbative scheme since Aju^ will 
turn out to be not only much smaller than the Newtonian monopole term, 
but also of the Schwarzschild-like planetary accelerations which are of the 
order of 

Agtr ^ = 1 X 10-9 - 6 X 10-11 m s-^. (6) 

c r 

Anyway, it would be, in practice, useless since the only addition with respect 
to the orbital effects resulting from the standard perturbative scenario would 
consist of further, small cross GTR-RIN orbital effects which, in the present 
case, would be of the order of Aju^^v / c)"^ . Since typical planetary speeds 
are of the order of ~ \jGMjr ~ 10^ m s~i, such mixed terms would be 
completely negligible. Let us remark that putting a bound smaller than the 
GTR Schwarzschild-type term is, actually, consistent. Indeed, GTR is fully 
modelled in the softwares routinely used to process astronomical data: thus, 
what actually is relevant is the uncertainty in the GTR terms, which is, at 
present, at a 0.01 percent level or less (see Section [3TT] below) . Moreover, as 
it will be shown in Section 13. 1|, the RIN and GTR effects have also different 
time signatures, thus allowing for tighter constraining of RIN. 



In the following we will neglect it to make the notation less cumbersome. 
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For analytical purposes, we will work in the R — T — N formalism by 
expressing the position perturbation Ar in terms of its radial, transverse and 
out-of-plane (or normal) projections AR, AT, AN onto the three orthogonal 
directions of the co-moving frame with unit vectors r, f , z), so that 

Ar = ARf + ATf + ANi'. (7) 

In order to analytically work out the range and range-rate perturbations 
due to some dynamical effects, let us review in some details some key- features 
of the R — T — N formalism. In terms of the standard Keplerian orbital ele- 
ments and of the usual unit vectors i, j, k of an inertial frame with Cartesian 
rectangular coordinates having its origin in the central body, the R — T — N 
versors, evaluated onto the unperturbed orbit, are [11] 

cos cos u — cos / sin sin u 
sin $7 cos u + cos / cos Q, sin u | (8) 
sin / sin u 

— sin u cos n — cos / sin cos u 

— sin sin u + cos / cos Q cos u | (9) 
sin / cos u 

sin / sin 

- sin/ cos I (10) 
cos /, 

where /, ri, uj are the inclination of the orbit to the reference {x, y} plane 
adopted, the longitude of the ascending node and the argument of the peri- 
centre, respectively; u = uj + f is the argument of latitude, in which / is the 
true anomaly reckoning the instantaneous position of the test particle along 
its Keplerian ellipse. Thus, after having analytically worked out the per- 
turbations of the Keplerian orbital elements, it will be possible to calculate 
AR, AT, AN according to Ref. [H] 



AR = (0 Aa-acos/Ae + ^g|AA1, 



AT = a sin / 



1 4- ^ 

^ a(l-e^) 



Ae + r(cosIAn + Auj)+ (^) Vl^AM 



AN = r(siniiA/ — cosusin/Af]), 

(11) 

in which a,e,A4 are the semimajor axis, the eccentricity and the mean 
anomaly, respectively. 
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The Gauss equations for the variation of the Keplerian orbital elements 
are [T3l 



dt 



de 
dt 



dJ 

dt 



dO, 
dt 



duj 
dt 



dM 
dt 



[eAnsmf + Ari^,)], 



na 



VT 



f{Ansmf + AT [cos/ + i(l-i)]}, 
^ An{^)cosu, 



1 



(12) 



na sin I\/l—e^ 



smu, 



ArcosJ + At + n sin f 



cos If, 



n 



In eq. (jl2p p = a(l — e^) is the semi-latus rectum, n = \J GMjo? is the 
unperturbed Keplerian mean motion related to the Keplerian orbital period 
by n = 27r/Pb) and Ar, At, Ai^f are the radial, transverse and out-of-plane 
components of the disturbing acceleration A which have to be computed 
onto the unperturbed Keplerian ellipse. It turns out that, in order to make 
the calculations easier, it is more convenient to use the eccentric anomaly E 
instead of the true anomaly /; basically, E can be regarded as a parametriza- 
tion of the polar angle in the orbital plane. To this aim, useful conversion 
relations are |13| 

cos E—e 



cos/ 
sin/ 
r 
dt 



1— e cos E ' 

Vl— e-^ sini? 
1— ecosi? ' 

o(l — e COS E), 
(hzecos_E\dE. 



(13) 



In the specific case of RIN, the R — T — N components of its acceleration 
are simply 

Ar = Aptin, 



At 
An 



0, 
0, 



(14) 
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in which ^Rin can be either positive or negative. Inserting eq. (jl4p into eq. 
(jl2p and integrating it by means of eq. ()13p from the initial value of the 
eccentric anomaly Eq to a subsequent, generic value E yield 



f Aa 

Ae 
AI 

An 

cos I An + Auj 
AM 



2eAain(cos _E— cos Eq) 



^Rin ( 1 ~ ) (cos — COS _Eo ) 



— £"0) — (^^^^~^^^^q) 



T" 



_3 _ + (£i£:E-sinEo) ^ 



e^(sin2E-sm2go) 



(15) 



+ 3e (sinS — sinE'o) 

Thus, eq. (fT5|) inserted into eq. pT]) yield the R — T — N position 
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perturbations due to RIN 
/ AE> ^Rin r 

~ n^(l-ecos£) t 

[cos {E - Eo) - 1] + 
+ 3e [{E - Eo) sin E + cosE - cos Eq] + 
+ Se^ [cos(£^ - £;o) - 1] + 

+ ^ [3 cos S + cos 3E -4 cos Eq (cos 2£; + sin E sin £'o)] 

AT^ ^Rin r 

~ n2(l-ecos£) t 

2 [sin {E - Eo) -{E- Eo)] + 
+ § [4 - £;o) cos £; - 6 sin ^ + sin {2E - Eq) + 5 sin Eq] - 
- [{E - Eo) cos^ E - (sin S - sin Eo) cos E^o] 
}, 

AN = 0. 

(16) 

Note that the results of eq. ([15]) and of eq. ([16]) are exact; no approximations 
in e have been used. Moreover, eq. (jlOp does not present any singularities 
for particular values of e. 

In order to conveniently plot eq. ()16|) as a function of time we will use 
suitable partial sums of the serie^ 

oo 2 



E = M + y^ -Jq{qe) sin (qM), (17) 
where Jq{qe) is the Bessel function of the first kind [15]. Indeed, the mean 



*It converges for all e < 1 like a geometric series with 
ratio r = [e/ (l + Vl - e^)] exp (Vl - e ^) [H]. See also 



http://mathworld.woliiam.com/KeplersEquation.htmH 
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anomaly is a parametrization of time according to 



M = n{t-tp), (18) 

where tp is the time of the passage at pericenter. It turns out that for 
small eccentricities just a few terms in eq. p7|) have to be retained. Note 
that eq. (fT6]) tells us that, for a given value of ^Rin, the largest effects 
occur for those orbiters having the largest orbital periods; in particular, for 
a given central body of mass M, the most distant test particles orbiting 
it experience the largest perturbations. Indeed, the shifts of eq. (jl6p are 
proportional to A^i^iPbi expected. This fact is important because, for 
a given level of accuracy in determining the orbits of the probes used, the 
tightest constraints on vlRin come just from the most distant ones with 
respect to M. 

The perturbations of eq. (fT6]l can fruitfully be used to analytically work 
out the two-body range and range-rate perturbations between two test parti- 
cles A and B orbiting the same central body of mass M. Indeed, concerning 
the range, from 

yo^ = irA-rB)-{fA-rB), 

(19) 

P p ' 

to be evaluated onto the unperturbed Keplerian ellipses of the two test 
particles A and B, it follows that, for a generic perturbation, the range shift 
Ap is [IE] 

Ap = (AfA - Are) • p. (20) 
In the case of the unperturbed Keplerian ellipse it is 

f = r f, (21) 

with r as in eq. (fT3]l . and r given by eq. ([8]). 

2.2 The two-body range-rate 

It is also possible to analytically work out the two-body range-rate perturba- 
tion Ap |16] by, first, working out the R — T — N shifts of the test-particle's 
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velocity. In general, they are [12 



AVR = (^eAa ^ ofAej - (cos /Al] + A^) - ^ A.M , 

AvT = -iHv^Aa + ^^jg^Ae + ^^J|^(cos/AO + Aw), 
Avn = ^""gi [(cos u + e cos a;) AI + (sin + e sin u) sin /All] 



so that 



(22) 

A-f; = Ai)R f + Avt t + jv z^- (23) 
In the case of RIN, inserting eq. (|15p in eq. (|22p straightforwardly yields 

^'^^ ~ n(l-ecos£;)2 t 

(/; - £;o) (2 + ecos E) - sin (/; - Eq) - 
— e {4 — ecos E) {s'mE — s'mEo) + 

+ ^ [e cos sin 2/; + (1 - e cos E) {2E - 2Eq + sin 2E) - sin 2£;o] 



Avt = [ 

^ n{l—ecosE) 

{E - Eq) sin E-1 + cos{E - Eo)+ 
+ e (cos — cos -Eq) 



t AvN = 0. 

(24) 

Also eq. (j24p are exact in e. Note that, for a given value of ^Rin, also 
the shifts of eq. (j24p get larger for more distant, i.e. slower, orbiting test 
particles around a central body; indeed, as expected, they are proportional 

to ^Rin-Pb- 

In order to work out the two-body range-rate perturbation, the following 
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unit vector, computed onto the unperturbed Keplerian ellipse, is needed [16] 



Pn = , (25) 

P 

where 

p={vA- vb) ■ p. (26) 

In eq. (|25p -eq. ^ thcl Keplerian test particle's velocity must be used; it 
is 

V = vn r + Vt t, (27) 

with 



VR 



Vt 



nae sin / 

(28) 

na(l+e cos /) 



Note that, by construction, pn is orthogonal to p. The two-body range-rate 
perturbation Ap is, thus, [16] 

Ap = [AvA - Avb) ■ p + (AfA - Are) • pn- (29) 

In the present specific case, inserting eq. (fT6]l and eg. (pl|) into eq. (f29]l . with 
eq. pT|) . allows to plot the range-rate perturbation due to RIN as a function 
of time. Finally, let us note that an alternative approach to compute Ap 
consists of straightforwardly taking the derivative Ap with respect to t after 
that its time series has been generated from eq. (j20p with the aid of eq. 



3 Confrontation with the observations 



3.1 The range and range-rate residuals 

Here we plot the analytically computed time-series of the two-body range 
and range-rate perturbations caused by RIN for A=Earth and B given by 
various bodies of the solar system. Then, we compare them with the existing 
residuals produced with the latest ephemerides in which RIN was not mod- 
eled in order to preliminarily put constraints on the magnitude of ^Rin- In 

^Here we drop the subscript U denoting the unperturbed trajectory. 
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order to unambiguously and correctly interpret our results, it should be no- 
ticed that all standard Newtonian and Einsteinian gravitational effects have 
been accurately modeled in computing the range-residuals, so that they 
account, in principle, for any unmodeled dynamical effect like just RIN. In 
particular, the modeling of GTR dynamical effects has reached a sub-percent 
level of accuracy. Indeed, reasoning in terms of the PPN parameters /3 and 
7, latest determinations of them from fitting dynamical models of different 
ephemerides (INPOP, EPM, DE) to large data sets covering about one cen- 
tury by independent teams of astronomers are accurate at ~ 10~^ — 10~^ 
level [THl [ig [201 [211 [231 [221 [Ml Ea [26]. We are confident in our results 
because we independently checked our analytically-produced time series by 
numerically integrating the equations of motion for the pairs A-B with and 
without RIN; the resulting range and range-rate signatures coincide with 
those analytically computed. 

According to Table 1 of Ref. |19j . the Mercury range residuals cover 
rather continuously a time span 33 yr long (1964-1997) with a Root-Mean- 
Square (RMS) error of 575 m; see also Figure B-2 and Figure B-3 of Ref. [28] . 
Figure [1] shows that the largest admissible value for ^Rin yielding an unmod- 
elled range signal still compatible with the existing residuals is 10~^^ m s~^; 
indeed, a larger value for A^nn would induce a signature larger than the 
actual residuals. 



Earth-Mercury range 




5 10 15 20 25 

1 (yr) 



Figure 1: Earth-Mercury range perturbation Ap, in m, caused by ^Rin = 

1 X 10~^^ m s^^ over At = 33 yr. The initial conditions have been retrieved 
from the NASA-JPL WEB interface HORIZONS. 

Table 1 of Ref. [19j tells us that the Venus range residuals span rather 
continuously an interval 34 yr long (1961-1995) with a RMS error of 584 m; 
see also Figure B-6 of Ref. |28j . The Magellan range-rate residuals covering 

2 yr (1992-1994) have a RMS error of just 0.007 mm s~^ [l9]. According 
to Figure [21 also in this case we have I^Rml ^ 10~^^ m s~^. In Figure [3| 
we depict the nominal Earth- Venus range-rate perturbation due GTR: cfr. 
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Earth-Venus range Eai'th-Venus range-rate 




Figure 2: Left panel: Earth- Venus range perturbation Ap, in m, caused 
by Ayhu = 2 X 10~^^ m s~^ over At = 34 yr. Right paneh Earth- Venus 
(Magellan) range-rate perturbation Ap, in mm s~^, caused by A^nn = 2 x 
j^Q-13 g-2 Q^gj, At = 2 yr. The initial conditions have been retrieved from 
the NASA-JPL WEB interface HORIZONS. 



with the right panel of Figure [2j As anticipated in Section I2.H the different 
patterns are evident. 



Earh- Venus range-rate 




Kyrl 



Figure 3: Earth- Venus range-rate nominal perturbation Ap, in mm s , 
caused by GTR over At = 2 yr. The initial conditions have been retrieved 
from the NASA-JPL WEB interface HORIZONS. 

Tighter constraints come from Mars. According to Table 1 of Ref. [19] . 
the 6-yr long range residuals of the Odyssey spacecraft (2002-2008) have a 
rms error of 1.2 m; see also Figure B-11 of Ref. [20]. Figure [5] shows that, in 
this case, ^Rin is constrained at a 10~^^ m s~^ level. Similar results come 
from the range-residuals of the Mars Global Surveyor (MGS) spacecraft, 
covering 8 yr (1998-2006) and accurate to 1.4 m [19]; Figure [4] shows that 
they practically constrain ^Rin at the same level. It turns out that the 
range-rate residuals of Viking and Pathfinder [19] do not yield constraints 
competitive with those from the range. The pattern of the GTR signal is 
different with respect to the RIN one, as clearly shown by Fig. 23 of Ref. 
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Earth-Mars range Ear lh-Mar s range 
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Figure 4: Left panel: Earth-Mars (Odyssey) range perturbation Ap, in m, 
caused by A^nn = 1 x 10"^^ m s~^ over At = 6 yr. Right panel: Earth-Mars 
(MGS) range perturbation Ap, in m, caused by A^nn = 1 x 10^^^ m s~^ over 
At = 8 yr. The initial conditions have been retrieved from the NASA-JPL 
WEB interface HORIZONS. 



|27| covering 5 yr. 

The present-day available range-residuals of the Cassini spacecraft or- 
biting Saturn cover about 2 yr (2004-2006) and are accurate to 20 m; see 
Table 1 of Ref. [19] and Figure B-20 of Ref. [28]. They allow to constrain 
^Rin at a 10~^^ m s~^ level, as shown by Figure [5l Also in this case, the 



Earth-Saturn range 




t(yr) 

Figure 5: Earth-Saturn (Cassini) range perturbation Ap, in m, caused by 
^Rin = 1 X 10~^^ m s~^ over At = 2 yr. The initial conditions have been 
retrieved from the NASA-JPL WEB interface HORIZONS. 

Earth-Saturn range signal due to GTR is quite different with respect to the 
RIN one, as it can be noticed from Fig. 45 of Ref. [27] covering 5 yr. 

At this point it is important to stress that our results clearly show that 
^Rin cannot be the cause of the Pioneer anomaly. Indeed, Arih is a con- 
stant for a specific system in consideration [1] in the sense that for a given 
central body of mass M acting as source of the gravitational field ^Rin is 
fixed: neither spatial nor temporal variations are admitted. Thus, it is not 
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possible that ^Rin, being constrained to 10~ — 10^ m s^" level by the 
inner planets, suddenly jumps to 10"^*' m s~^ in the regions in which the 
Pioneer anomaly manifested itself in its presently known form. See also the 
discussion at the end of Section 13.21 in which different empirical quantities 
are used to constrain Ap^o after 20 au. 

Moving to the neighborhood of the Earth, Figure B-1 of Ref. [28] shows 
that the residuals of the Earth-Moon range, constructed from the data con- 
tinuously collected with the Lunar Laser Ranging (LLR) technique [29], are 
at a cm-level since about 1990. Figure [6] depicts the lunar range signature 
over At = 20 yr caused by a terrestrial RIN ^Rin = 5 x 10~^^ m s~^; its 
magnitude is as large as 4 cm at most, in agreement with the lunar laser 
ranging residuals. Larger values for ^Rin would yield signals too large; for 
example, it turns out that Aj^n = 1 x 10~^^ m s~^ yields a RIN range signal 
of almost 10 cm. 



Earth-Moon range 




1000 2000 3000 4000 5000 6000 7000 



Figure 6: Earth- Moon range perturbation Ap, in cm, caused by a terrestrial 
^Rin = 5 X 10~^^ m s~^ over At = 20 yr. The initial conditions have been 
retrieved from the NASA-JPL WEB interface HORIZONS. 

It must be pointed out that, in principle, the bounds obtained so far 
might be somewhat too tight since the unmodeled RIN signatures, if they 
really existed in the data, may have been partially removed in the estimation 
of the initial conditions in the data reduction process; actually, RIN should 
be explicitly modeled, and a dedicated parameter accounting for it should be 
solved-for along with the other ones routinely estimated in the usual way. 
First attempts towards the implementation of such a strategy applied to 
a generic extra-acceleration A radially directed towards the Sun may have 
been recently performed by Folkner in Ref. |20j . He obtained from the DE 
ephemerides by JPL an upper bound of A < 1 x 10~^^ m s~^ on the Earth- 
Cassini range [20], in agreement with our Figure [5j On the other hand, 
the bound obtained from Mars in Ref. [20] is A < 8 x 10~^^ m s~^, which 
is about one order of magnitude larger than ours in Figure 21 However, 
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it must be noted that, basically, no details have been released in Ref. [20] 
either concerning the data set of the spacecrafts used for Mars nor about 
the methodology adopted to constrain A from them. 

3.2 The perihelion precessions 

Another possible approach to constrain A^i^i consists of looking at the em- 
pirically determined corrections Aw to the standard Newtonian/Einsteinian 
secular precessions of the longitudes of the planetary perihelia zu = Q + ui, 
and compare them to theoretically predicted RIN-type precessions of w |17] . 

The corrections Azu have recently been estimated as solved-for parame- 
ters by independent teams of astronomers [HI [211 [22] by fitting accurate 
dynamical force models of their ephemerides, which only include the usual 
Newtonian/Einsteinian dynamics, to long observational records spanning al- 
most one century. They are listed in Table [TJ 

The secular precession of w caused by a constant and uniform extra-acceleration 
A radially directed towards the Sun has been computed by the authors 
of Refs. |30[ [3T| [32] in the framework of the investigations of the Pioneer 
anomaly. It is 

^ = -A^^. (30) 
dt na 



By applying eq. (j30p to RIN, a comparison of it with the most accurate 
values of Avu in Table [1] yields the upper bounds on Aju^ listed in Table [2j 
Remarkably, the constraints on ^Rin of Table [21 obtained from the plane- 
tary perihelia, are in substantial agreement with those of Section [3TT] inferred 
from the range and range— rate residuals. The perihelion of the Earth yields 
a bound as low as ^Rin < 7 x 10~^^ m s~^. 

Incidentally, it can also be noticed from Table [2] that the largest admis- 
sible values for ^Rin from the perihelia of Uranus, Neptune and Pluto are 
smaller than the smallest value of the anomalous Pioneer acceleration, i.e. 
^Pio = 7.41 X 10~^° m s~^ (IcT-level), by a factor 18 - 25, respectively. At 
3(T— level, i.e. for Apio = 4.75 x lO^^'' m s^^, they are smaller than it by a 
factor 12 - 16. 
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Table 1: Estimated corrections Act, in milliarcseconds per century (mas 
cty~^), to the standard Newtonian/Einsteinian secular precessions of the 
longitudes of the perihelia ct of the eight planets plus Pluto determined 
with the EPM2008 [H], the INPOP08 [21], and the INPOPlOa [22] 
ephemerides. Only the usual Newtonian/Einsteinian dynamics was mod- 
elled, so that, in principle, the corrections Act account for any other un- 
modelled/mismodelled dynamical effect. Concerning the values quoted in 
the third column from the left, they correspond to the smallest uncertainties 
reported in Ref. |21j . Note the small uncertainty in the correction to the 
precession of the terrestrial perihelion, obtained by processing Jupiter VLBI 
data [21]. 



Planet 


Act [ig 


Act t2lj 


Act [22] 


Mercury 


-4 ±5 


-10 ±30 


0.2 ±3 


Venus 


24 ±33 


-4 ±6 




Earth 


6±7 


0± 0.016 




Mars 


-7±7 


0±0.2 




Jupiter 


67 ±93 


142 ± 156 




Saturn 


-10± 15 


-10 ±8 


0±2 


Uranus 


-3890 ± 3900 


± 20000 




Neptune 


-4440 ± 5400 


± 20000 




Pluto 


2840 ± 4510 







4 Summary and conclusions 

We analytically worked out the perturbing effects which a Rindler-type 
anomalous acceleration ^Rin, naturally arising from a 2-dimensional dilaton- 
based long-range modification of gravity, would induce on the orbital motion 
of a test particle orbiting a central body acting as source of the modified 
gravitational field. In particular, in view of a comparison with the most re- 
cent observations we focussed on the effects of Aju^ on the two-body range 
p and range-rate p because they are direct, unambiguous observables widely 
used in satellite and planetary investigations. 

It turns out that ^Rin does actually affect p and p with long-term signa- 
tures which can fruitfully be compared with the residuals for such observ- 
ables built by processing extended data records for some planets and the 
Moon with the latest ephemerides in which the Rindler perturbation has 
not been explicitly modeled. It turns out that Mercury and Venus constrain 
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Table 2: Upper bounds on the magnitude of RIN, in m s^^, from a com- 
parison of the theoretical prediction of eq. ()30p for the RIN-induced secular 
precession of the longitude of the pericenter w of a test particle and the most 
accurate values of the empirically determined corrections Aro to the stan- 
dard Newtonian/Einsteinian secular planetary precessions of the perihelia 
quoted in Table [TJ 



Planet 


^Rin ( 


m g 




Mercury 


< 2 X 


10" 


13 


Venus 


< 3 X 


10-- 


-13 


Earth 


< 7 X 


10- 


-16 


Mars 


< 7 X 


10- 


-15 


Jupiter 


< 2 X 


10" 


-12 


Saturn 


< 3 X 


10- 


-14 


Uranus 


< 4 X 


10- 


-11 


Neptune 


< 4 X 


10- 


11 


Pluto 


< 3 X 


10- 


11 



a solar ^Rin to a level of about IQ-^'^ m s"^. The bounds inferred from recent 
data sets of some spacecrafts orbiting Mars are of the order of 10"^^ m s~^, 
while the Cassini range residuals yield ^Rin ^ 1 x IQ-^^ m s"^ at the Sat- 
urn's distance. Analogous constraints are obtained by using the empirically 
determined corrections A-cu to the standard Newtonian/Einsteinian secular 
precessions of the planetary perihelia: the Earth yields a tighter bound as 
low as ^Rin ^ 7 X IQ-^^ m s~^. The tightest constraints come from the lunar 
range residuals according to which ^Rin ^ 5 x 10~^^ m s"^ for a modifica- 
tion of the Earth's gravitational field. Another analysis existing in literature, 
likely based on a different method, points towards ^Rin < 8 x IQ-^^ m s~^ 
at the Mars' distance, while it agrees with our result for Saturn. 

In principle, ^Rin should be explicitly modeled, the entire planetary and 
satellite data sets used should be re-processed with such modified dynamical 
models, and a dedicated solve-for parameter should be estimated. Such an 
approach is certainly rather expensive in terms of computational burden 
and time required, especially if other, more conventional tasks are pressing 
and have to be mandatorily performed. Moreover, if one is interested in 
other exotic effects predicted by some different modified gravities the entire 
procedure has to be repeated with the new model. Instead, the approach 
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followed here in the case of the Rindler-type modification of gravity at large 
distances can be easily and straightforwardly extended to other long-range 
modified models of gravity as well 
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